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We search for new charmless decays of neutral 6-hadrons to pairs of charged hadrons with the 
upgraded Collider Detector at the Fermilab Tevatron. Using a data sample corresponding to 1 fb _1 
of integrated luminosity, we report the first observation of the _B° — ► K~-k + decay, with a significance 
of 8.2a, and measure B(B° -> K~n + ) = (5.0 ± 0.7 (stat) ± 0.8 (syst)) x 10~ 6 . We also report the 
first observation of charmless b-baryon decays in the channels — > pir~ and A" — ► pK~ with 
significances of 6.0a and 11.5(7 respectively, and we measure B(A° — > piv~) = (3.5 ± 0.6 (stat) ± 
0.9 (syst)) x 10~ 6 and B(Ag -» pK~) = (5.6 ± 0.8 (stat) ± 1.5 (syst)) x 10~ 6 . No evidence is 
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7r + 7r~) < 1.2 x 10~ 6 at the 90% confidence level. All quoted branching fractions are measured using 
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Two-body non-leptonic charmless decays of 6-hadrons 
are among the most widely studied processes in flavor 
physics. The variety of open channels involving similar 
final states provides crucial experimental information to 
improve the accuracy of effective models of strong inter- 
action dynamics. The quark-level transition b — > u makes 
decay amplitudes sensitive to 7, the least known an- 
gle of the quark-mixing (Cabibbo-Kobayashi-Maskawa, 
CKM) matrix. Significant contributions from higher- 
order ('penguin') transitions provide sensitivity to the 
possible presence of new physics in internal loops, if the 
observed decay rates are inconsistent with expectations. 

Rich experimental data are currently available for B + 
and B° mesons, produced in large quantities in T(4S') de- 
cays [![, while much less is experimentally known about 
the charmless decay modes of the B°, which are expected 
to exhibit an equally rich phenomenology Information 
from B® decays is needed to better constrain the phe- 
nomenological models of hadronic amplitudes in heavy 
flavor decays. This would lead to increased precision 
in comparing data to predictions, allowing extraction of 
CKM parameters from non-tree-level amplitudes [2J and 
greater sensitivity to new physics contributions. 

Of the possible B® decay modes into pairs of charmless 
pseudoscalar mesons, only the B® — > K + K~ has been 
observed to date The B® — > K~tt + is of particular 
interest, because its branching fraction is sensitive to the 
CKM angle 7 0] and the current experimental bound Q 
is lower than most predictions [!, H, @] . 

A measurement of the branching fraction of the B® — > 
7r + 7r~ mode, along with the B° — » K + K~ mode, 
would allow a determination of the strength of penguin- 
annihilation amplitudes which is currently poorly 
known and a source of significant uncertainty in many 
calculations @. The present search is sensitive to both 
modes. Two-body charmless decays are also expected 
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from bottom baryons. The modes A° — > pK~ and 
A° — > pn~ are predicted to have measurable branching 
fractions, of order 10~ 6 9], and, in addition to the inter- 
est in their observation, must be considered as a possible 
background to the rare B° B and B° modes being investi- 
gated. 

In this Letter we report the results of a search for rare 
decays of neutral bottom hadrons into a pair of charged 
charmless hadrons (p, K or 7r), performed in 1 lb -1 of pp 
collisions at y/s = 1.96 TeV, collected by the upgraded 
Collider Detector (CDF II) at the Fermilab Tevatron. 
We report the first observation of modes B® — > K~tt + , 
A® — > pK~ , and A° — * pir~ , and measure their relative 
branching fractions (Toj . 

CDF II is a multipurpose magnetic spectrometer sur- 
rounded by calorimeters and muon detectors. The detec- 
tor components relevant for this analysis are briefly out- 
lined below; a more detailed description can be found in 
Ref. A silicon microstrip vertex detector (SVX) and 
a cylindrical drift chamber (COT) immersed in a 1.4 T 
axial magnetic field allow reconstruction of charged- 
particle trajectories (tracks) in the pseudorapidity range 
1 77 1 < 1.0 [12j. The SVX consists of six concentric lay- 
ers of double-sided silicon sensors with radii between 
2.5 and 22 cm, each providing a measurement with up 
to 15 (70) /Ltm resolution in the <f> (z) direction. The 
COT has 96 measurement layers, between 40 and 137 
cm in radius, organized into alternating axial and ±2° 
stereo superlayers. The transverse momentum resolution 
is a PT jp\ ~ 0.15%/(GeV/c), corresponding to a typical 
mass resolution of 22 MeV/c 2 for our signals. The specific 
ionization energy loss (dE / dx) of charged particles in the 
COT can be measured from the collected charge, which 
is logarithmically encoded in the output pulse width of 
each wire, and provides 1.5a separation between kaons 
and pions with momenta greater than 2 GeV/c. 

The data were collected by a three-level trigger sys- 
tem, using a set of requirements specifically aimed at se- 
lecting two-pronged B decays. At level 1, COT tracks 
are reconstructed in the transverse plane by a hard- 
ware processor (XFT) [13]. Two opposite-charge par- 
ticles are required, with reconstructed transverse mo- 
menta Pti,Pt2 > 2 GeV/c, the scalar sum p?i + Pt2 > 
5.5 GeV/c, and an azimuthal opening-angle A(j) < 135°. 
At level 2, the silicon vertex trigger (SVT) [lj| combines 
XFT tracks with SVX hits to measure the impact pa- 
rameter d (distance of closest approach to the beam line) 
of each track with 45 /jm resolution. The requirement of 
two tracks with 0.1 < d < 1.0 mm reduces the light quark 
background by two orders of magnitude while preserving 
about half of the signal. A tighter opening-angle require- 
ment, 20° < Acf> < 135°, preferentially selects two-body 
B decays over multi-body decays with 97% efficiency and 
further reduces background. Each track pair is then used 
to form a B candidate, which is required to have an im- 
pact parameter ds < 140 /im and to have travelled a dis- 
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tance Lt > 200 /im in the transverse plane. At level 3, 
an array of computers confirms the selection with a full 
event reconstruction. The overall acceptance of the trig- 
ger selection is ~ 2% for 6-hadrons with pt > 4 GeV/c 
and \r]\ < 1. 

The offline selection is based on a more accurate de- 
termination of the same quantities used in the trigger, 
with the addition of two further observables: the iso- 
lation (Is) of the B candidate (l5(, and the quality of 
the three-dimensional fit (x 2 with 1 d.o.f.) of the decay 
vertex of the B candidate. Requiring a large value of 
1b reduces the background from light-quark jets, and a 
low x 2 reduces the background from decays of different 
long-lived particles within the event, owing to the good 
resolution of the SVX detector in the z direction. The 
selection is optimized for detection of the B® — > K~w + 
mode. Maximal sensitivity for both discovery and limit 
setting is achieved with a single choice of selection re- 
quirements [l6| by minimizing the variance of the esti- 
mate of the branching fraction in the absence of signal 
171 ]. The variance is evaluated by performing the full 
measurement procedure on simulated samples contain- 
ing background and all signals from the known modes, 
but no B® — ► K~ir + signal. The background fraction for 
each selection is determined from data by extrapolating 
the mass sidebands of the signal, and the signal yield is 
predicted by a detailed detector simulation. This pro- 
cedure yields the final selection: 1b > 0.525, x 2 < 5, 
d > f 20 fim, ds < 60 fim, and Lt > 350 fim. 

No more than one B candidate per event is found af- 
ter this selection, and a mass (m W7r ) is assigned to each, 
using a charged pion mass assignment for both decay 
products. The resulting mass distribution is shown in 
Fig. [U A large peak is visible, dominated by the over- 
lapping contributions of the B° — ^ K + ir~ , B a — > 7r + 7r~ , 
and B a s K+R- modes 0. A B° K+R- signal 
would appear as an enhancement around 5.f8 GeV/c 2 , 
while signals for the other modes of this search are ex- 
pected at masses higher than the main peak (5.33-5.55 
GeV/c 2 ). Backgrounds include mis-reconstructed multi- 
body 6-hadron decays (physics background) and random 
pairs of charged particles (combinatorial background) . 

We used an unbinncd likelihood fit, incorporating kine- 
matic (kin) and particle identification (PID) information, 
to determine the fraction of each individual mode in our 
sample. The likelihood for the ith event is 



d = (l-&)^/,£f\C 



PID 

3 



kin /*PID\ 



b{f p £^ n £ p p lu + (l-f p )£f n £ 



(1) 



where the index j runs over all signal modes, and the 
index 'p' ('c') labels the physics (combinatorial) back- 
ground terms. The fj are the signal fractions to be deter- 
mined by the fit, together with the background fraction 
parameters b and f p . 




5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 

Invariant rat-mass [GeV/c 2 ] 

FIG. 1: Mass distribution of reconstructed candidates. The 
charged pion mass is assigned to both tracks. The total pro- 
jection and projections of each signal and background compo- 
nent of the likelihood fit are overlaid on the data distribution. 
Signals and multi-body B background components are shown 
stacked on the combinatorial background component. 



The kinematic information is summarized by three 
loosely correlated observables: (a) the mass rn W7V ; (b) the 
signed momentum imbalance a = (1— pi/p2)qi, where pi 
(P2) is the lower (higher) of the particle momenta, and q\ 
is the sign of the charge of the particle of momentum p\ ; 
(c) the scalar sum of particle momenta ptot — Pi+P2- The 
above variables allow evaluation of the invariant mass 
77ii2 of a candidate for any mass assignment of the decay 
products (777,1,777,2), using the equation 

m i2 = m lir - 2m l + m\ + m 2 2 + 
-1\jp\ + ml\jp 2 2 + ml + 2^p 2 + m\^jp\ + m 2 , (2) 



where pi = ijzj^Ptot , Pi 
We used 
[5.00,5.12] U [5.6,6.2] GeV/c 2 



iPtot- 



2-|a| 

sidebands in data (m^ x e 
z 2 ) to obtain the kinematic 
distributions of backgrounds [13] • The mass distribu- 
tion of the combinatorial background is parameterized 
by an exponential function, while the physics background 
is modeled by an ARGUS function I8( convoluted with 
a Gaussian resolution function. In order to ensure the 
reliability of the search for small signals in the vicin- 
ity of larger peaks, the shapes of the mass distribu- 
tions assigned to each signal have been modeled in de- 
tail. We have included the momentum dependence and 
non-Gaussian tails of resolution from a full simulation 
of the detector, and the effects of soft photon radiation 
in the final state, based on recent QED calculations [lj| . 
This resolution model was checked against the observed 
shape of the D° — * K~ir + signal in a sample of 1.5 x 10 6 
D* + — > _D°7r + decays, collected with a similar trigger se- 
lection. The observed discrepancies are below the I0 -3 
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FIG. 2: Distribution of the relative signal likelihood, Cs/(£s + £other), in the region 5.1 < < 5.6 GeV/c 2 . For each event, 
C s is the likelihood for the B° -» K~tt + (a), A° b -> pK" (b), or A° b — > p7T (c) signal hypotheses, and £ ther is the likelihood for 
everything but the chosen signal, i.e. the weighted combination of all other components according to their measured fractions. 
Points with error bars show the distributions of data and histograms show the distributions predicted from the measured 
fractions. 



level, and their effect on the present measurement is neg- 
ligible in comparison with other systematic uncertainties. 
The D*+ -> D°tt + sample was also used to calibrate the 
dE I dx response of the drift chamber to kaons and pions, 
using the charge of the D* + pion to identify the D° de- 
cay products. The dE/ dx response of protons was deter- 
mined from a sample of about 124,000 A — > pTr~ decays. 
The model of the background allows for pion, kaon, pro- 
ton, and electron components, whose fractions are deter- 
mined by the fit. Muons are indistinguishable from pions 
with the available 10% fractional dE / dx resolution and 
are therefore incorporated into the pion component. 

From the signal fractions returned by the likelihood fit 
we calculate the signal yields shown in Table U The sig- 
nificance of each signal is evaluated as the ratio of the 
yield observed in data, and its total uncertainty (statis- 
tical and systematic) as determined from a simulation 
where the size of that signal is set to zero. This evalua- 
tion assumes a Gaussian distribution of yield estimates, 
supported by the results obtained from repeated fits to 
simulated samples. This procedure yields a more accu- 
rate measure of significance with respect to the purely 
statistical estimate obtained from y/— 2Aln(£). We ob- 
tain significant signals for the B® — > K~ir + mode (8.2er), 
and for the A° -> pir~ (6.0a) and A° -> pK~ (11.5a) 
modes. Figure [2] shows relative likelihood distributions 
for these modes. No evidence is found for the modes 
— > 7r + 7r~ or B° — > K + K~ , in agreement with expec- 
tations of significantly smaller branching fractions. 

To avoid large uncertainties associated with produc- 
tion cross sections and absolute reconstruction efficiency, 
we measure all branching fractions relative to the B° — > 
K + ir~ mode. Frequentist upper limits [13] at the 90% 
C.L. are quoted for the unseen modes. For the mea- 
surement of A° branching fractions, the additional re- 



TABLE I: Yields and significances of rare mode signals. The 
first quoted uncertainty is statistical, the second is systematic. 



Mode 




N s Si f 


mificance 


fl" -» 


K~TV + 


230 ± 34 ± 16 


8.2a 




TT + TV" 


26 ± 16 ± 14 


< 3a 


B° —> 


K+K~ 


61 ± 25 ± 35 


< 3(7 






156 ± 20 ± 11 


11.5(7 




pn~ 


110 ± 18 ± 16 


6.0(7 



quirement pr(A%) > 6 GeV/c was applied to allow 
easy comparison with other A° measurements at the 
Tevatron, which are only available above this thresh- 
old [2^, [2^ . This additional requirement lowers the Ag 
yields by about 20%. The raw fractions returned by the 
fit were corrected for the differences in selection efficien- 
cies between different modes, which range from 8% to 
40% for the measurements of 6-mesons and A° branch- 
ing fractions, respectively. These corrections were deter- 
mined from detailed detector simulation, with the fol- 
lowing exceptions that were measured from data: the 
momentum-averaged relative isolation efficiency between 
B° and B°, 1.00 ± 0.03, has been determined from fully- 
reconstructed samples of B®— > J/tp<p, and B°^ J/ijjK*° 
decays 17]; the difference in efficiency for triggering on 
kaons and pions due to the different specific ionization in 
the COT (a « 5% effect) was measured from a sample 
of D + — ► K~n + Tt + decays triggered on two tracks, using 
the unbiased third track [24[ . Possible differences in effi- 
ciency of the isolation requirement between B° and Ag, 
and in the trigger efficiency between kaons and protons, 
were taken into account in the systematic uncertainties. 

The dominant contributions to the systematic uncer- 
tainty are the uncertainty on the combinatorial back- 
ground model and the uncertainty on the dE / dx cali- 
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TABLE II: Measured relative branching fractions of rare modes. The ratio /a/ fd is pr-dependent [2(j, and is defined here 
as: fh/fd = c(pp — > A°X;pr > 6 GeV/c, |^| < l)/a(pp — > B°X;pt > 6 GeV/c, |tj| < 1). Absolute branching fractions were 
derived by normalizing to the current world-average value B(B° -> K + tt~) = (19.4 ± 0.6) x 10~ 6 , and assuming the average 
values at high energy for the production fractions: f s /fd = 0.276 ± 0.034, and fk/fd = 0.230 ± 0.052 The first quoted 

uncertainty is statistical, the second is systematic. 

Mode Relative B Absolute £(10~°) 

B° -> K-TT+ h. ^^K~* l± Z 0.071 ± 0.010 ± 0.007 5.0 ± 0.7 ± 0.8 

B° -> TT+TT- & bIbSZk+C-) = - 007 ± - 004 ± - 005 - 49 ± - 28 ± °' 36 (< 12 at 90% C - L -) 

5° -> K + K" = °- 020 ± °- 008 ± °- 006 °- 39 ± °- 16 ± °- 12 (< °- 7 at 90% C - L 

AS -> pA"" S + »}3~*1+ = 0.066 ± 0.009 ± 0.008 5.6 ± 0.8 ± 1.5 

Ag -> P 7T- & b($SkI~J-) = °- 042 ± °- 007 ± °- 006 3.5 ± 0.6 ± 0.9 



bration and parameterization. Other contributions come 
from trigger efficiencies, physics background shape and 
kinematics, 6-hadron masses and lifetimes, and the pos- 
sible polarization of A° decays. 

The final results are listed in Table |TTJ Absolute 
branching fractions are also quoted, by normalizing to 
world-average values of production fractions and B(B° — > 
K + n~) [l|, |2l|. The branching fraction measured for the 
— > K~ir + mode is consistent with the previous upper 
limit (< 5.6 x 10~ 6 at 90% C.L.), based on a subsam- 
ple of the current data [3( • This agrees with the predic- 
tion in Ref. 1251 . but it is lower than most other predic- 
tions 0, |26| . The B® — > tt + ii~ upper limit improves 
and supersedes the previous best limit [3j. The present 
measurement of B(B Q — > K + K~) is in agreement with 
other existing measurements and has a similar resolu- 
tion [21] . but the resulting upper limit is weaker due 
to the observed central value. The sensitivity to both 
B° — > K + K~ and B® — > 7r + 7r _ is now close to the upper 
end of the theoretically expected range 0, 0, H3] • We 
also report the first branching fraction measurements of 
charmless A^ decays. They are significantly lower than 
the previous upper limit of 2.3 x 10 -5 [28|, and in rea- 
sonable agreement with predictions 0], thus excluding 
the possibility of large (O(10 2 )) enhancements from Im- 
parity violating supersymmetric scenarios [29J . Their ra- 
tio can be determined directly from our data with greater 
accuracy than the individual values. For this purpose, 
the additional px > 6 GeV/c requirement is not neces- 
sary, and we can exploit the full sample size, obtaining 
B(A° b -► pir-)/B(A° b -> pK~) = 0.66 ± 0.14 ± 0.08 in 
good agreement with the predicted range 0.60-0.62 [9(. 

In summary, we have searched for rare charmless de- 
cay modes of neutral 6-hadrons into pairs of charged 
hadrons in CDF data. We report the first observation 
of the modes B® — > K~tt + , A° — * pir~ , and A° — > pK~ , 
and measure their relative branching fractions. We set 
upper limits on the unobserved modes B° — ► K + K~ and 
B° s -> TT+TT- . 
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